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COMPUTER PROGRAM OF DATA REDUCTION PROCEDURES FOR 
FACILITIES USING C02-N2-02-Ar EQUILIBRIUM 
REAL-GAS MIXTURES 

By Charles G. Miller m 
Langley Research Center 

SUMMARY 

Data reduction procedures for determining free-stream and post-normal- shock flow 
conditions are presented. These procedures are applicable to flows of CO 2 , N 2 , O 2 , Ar, 
or mixtures of these gases and include the effects of dissociation and ionization. The 
assumption of thermochemical equilibrium free-stream and post-normal-shock flow is 
made. Although derived primarily to meet the immediate needs of the expansion tube of 
the Langley hot gas radiation research facility, these procedures are applicable to any 
supersonic or hypersonic test facility using these gases or mixtures thereof. 

The data reduction procedures are based on combinations of three of the following 
flow parameters measured in the immediate vicinity of the test section: 

(1) Stagnation pressure behind normal shock 

(2) Free-stream static pressure 

(3) Stagnation-point heat-transfer rate 

(4) Free-stream velocity 

(5) Free-stream density 

Thus, these procedures do not depend explicitly upon measured or calculated upstream 
flow parameters. The procedures are incorporated into a single computer program 
written in FORTRAN IV language. A listing of this computer program is presented, along 
with a description of the inputs required and a sample of the data printout. 

INTRODUCTION 

As discussed in reference 1, significant differences between measured test-section 
flow quantities and predicted test-section flow quantities based on upstream flow conditions 
were observed in the Langley pilot model expansion tube. In order to provide a means for 
obtaining accurate test-section conditions, a computational scheme based on flow proper- 
ties measured in the immediate vicinity of the test section was derived in reference 1. 



This scheme eliminated an explicit dependence upon measured or calculated upstream 
flow properties, thereby resulting in a substantial reduction in the uncertainty in pre- 
dicted test-section flow conditions. The study of reference 1 was, however, limited to 
thermochemical equilibrium, imperfect, real-air flows. 

Experimental investigations made in test facilities for which the test gas will be 
representative of Mars and Venus atmospheres (CC >2 predominance) require accurate 
predictions of test-section flow conditions for such real-gas mixtures. Therefore, the 
purpose of this program is to provide a means for obtaining accurate test- section condi- 
tions for real gases CO 2 , N 2 , 02, and Ar or mixtures of these gases. This is achieved 
by employing a computational scheme based on combinations of three flow parameters 
measured in the immediate vicinity of the test section. The measured flow properties 
serving as inputs to the data reduction procedures presented herein are as follows: 

(1) Stagnation pressure behind normal shock 

(2) Free-stream static pressure 

(3) Stagnation- point heat-transfer rate 

(4) Free-stream velocity 

(5) Free-stream density 

Of the 10 possible data reduction procedures from these inputs, seven are presented 
herein. (Data reduction procedures deemed less preferable in the error analysis study 
of ref. 1 are omitted.) The thermodynamic properties for CO 2 , N 2 , O 2 , Ar, and mixtures 
of these gases are obtained from references 2 and 3. Although the data reduction pro- 
cedures presented herein were derived to meet the immediate needs of the expansion tube 
of the Langley hot gas radiation research facility, these real -gas -mixture procedures are 
also applicable to supersonic or hypersonic test facilities using these gases or mixtures 
of these gases. 

The procedures are incorporated into a single computer program written in 
FORTRAN IV language. A listing of this computer program is presented in appendix A 
and a description of the inputs required is presented in appendix B. Program output 
parameters for various gas mixtures are illustrated by sample data printouts presented 
in appendix C. 


SYMBOLS 

The International System of Units (SI) is used for all physical quantities in this 
study. Conversion factors relating the International System of Units to the U.S. 
Customary Units are given in reference 4. 
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a 


speed of sound, m/sec 


h specific enthalpy, m2/sec2 (J/kg) 

M Mach number, U/a 

Np r Prandtl number 

N Re Reynolds number, m - * 

p pressure, N/m2 

q heat-transfer rate, W/m2 

R universal gas constant, 8.31434 x 10^ J/kmole-°K 


r nose radius, m 

sW 0/ /R nondimensional specific entropy 
T temperature, °K 

U velocity, m/sec 

W molecular weight, kg/kmole 

W 0 molecular weight of undissociated gas or gas mixture, kg/kmole 

X mole fraction, kmole of species i/kmole of mixture 


Yi 

Z 

1 


number of kmole of species i per mass of mixture, 
kmole of species i/kg of mixture 

compressibility factor, pW 0 /pRT 

isentropic exponent, 

parameter defined in equation (12), kg/'m3/2_ sec _ (N/m2)^^ 


’ajogjA 

3 !og P/sW 0 /R 
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viscosity, N-sec/m2 


M 

p density, kg/m 3 

Subscripts: 
c calculated 

g geometric 

i individual species 

m measured 

mix mixture 

prev previous value of a parameter 

t stagnation conditions behind normal shock 

w model wall 

1 free stream 

2 static conditions immediately behind normal shock 
Approximate value is denoted by 


ANALYSIS 

Before the discussion of the calculation procedures for determining free-stream 
and post-normal-shock flow quantities, the source of thermodynamic properties for 
C02-N2-02-Ar gas mixtures in chemical equilibrium employed in this study is discussed. 
Next, the iterative methods for crossing the standing normal shock are discussed, since 
these methods are common to all calculation procedures. Finally, the relations used to 
predict stagnation-point heat-transfer rate are discussed. 
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Thermodynamic Properties for Arbitrary Gas Mixtures 

Thermodynamic properties for C02-N2-02-Ar gas mixtures are obtained from the 
results of references 2 and 3. These results include dissociation and ionization. Basic 
assumptions used in obtaining these results are as follows: 

(1) The mixture is composed of ideal gases. 

(2) For diatomic species the rigid-rotor harmonic oscillator model is used with 
vibrational- rotational corrections for each electronic configuration. 

(3) Only electronic levels with principal quantum number less than or equal to five 
are included. 

The procedure of references 2 and 3 is based upon the free -energy minimization method 
of reference 5. For a given pressure and temperature, the free energy for individual 
species are computed from the partition function of statistical mechanics. (The five 
components and 26 species considered in refs. 2 and 3 are tabulated in appendix B.) The 
equilibrium composition is then obtained by minimization of the free energy. After the 
composition for a given pressure and temperature is determined, the corresponding ther- 
modynamic properties of interest are computed directly. 

The program of references 2 and 3 is included in the present study as a subroutine 
and referred to herein as ROGO. The basic inputs to ROGO, other than the atomic energy 
level constants (which are tabulated in ref. 2), are pressure, temperature, and initial 
estimates of the species concentrations. For the free-energy minimization method to 
determine the composition of a gas mixture, the initial estimates of Yj are somewhat 
arbitrary. In the present study, the partial pressures of the CO 2 , N 2 , O 2 , and/or Ar test 
gases in the initial state (prior to addition of energy) are assumed known. The corre- 
sponding mole fractions for these gases are obtained from Dalton's law. 


*i = 


Pi 

Pmix 


( 1 ) 


and the molecular weight is given by 

w o = 1 XjWj 
i 


( 2 ) 


The initial estimates for the species concentrations for CO 2 , N 2 , O 2 , and Ar are obtained 
from 


Yi = 


ii 

w 0 


(3) 
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For the remaining 22 species, the initial estimates for the concentrations are set equal 
to 10"20. Thermodynamic properties obtained from ROGO, for a given p and T, are 
a, h, sW 0 /r, Xi, Z, y E , and p. 

Combinations of input thermodynamic properties other than p and T are 
required in the present study. An iterative -interpolation scheme was derived to enable 
determination of thermodynamic properties from combinations of h, p, sW 0 /R, and p. 
This iterative -interpolation scheme is referred to herein as FIND (I), and the inputs for 
a given I are as follows: 


I 

Inputs to FIND (I) 

1 

P* P 

2 

P, h 

3 

p, sW 0 /R 

4 

h, p 


In the iterative phase of FIND (1, 2, or 3), T is varied (for a given p) to bring the cor- 
responding p, h, or sW 0 /r from ROGO toward agreement with the known p, h, or 
sW 0 /r. This iteration phase is terminated when the uncertainty in T is reduced to less 
than 250° K. ROGO is used to obtain the thermodynamic properties at the midpoint and 
ends of this region of uncertainty. An interpolation is performed over this region to 
obtain the desired thermodynamic properties. To assist in the reduction of computer 
time, the flow region (free stream or postshock) is identified in FIND (I). Upper and 
lower limits on T for these two regions are inputs. Hence, if the user of this program 
has a priori approximate values of Tj and T 2 or T^, a reduction in computer time 
can be realized by establishing these limits on T (see appendix B); if not, the limits are 
taken as 0° K to 8000° K for free-stream conditions and 0° K to 16 000° K for postshock 
conditions. 

For FIND (4), where h and p are inputs, an initial value of p is assumed and 
the iteration-interpolation scheme for FIND (1) employed to determine the thermodynamic 
properties for this p (that is, with this p and p, T is varied until the p from 
ROGO is within 0.1 percent of the input p). The calculated h is compared to the input 
h; if not within 0.1 percent, p is varied and the procedure repeated. This iterative 
scheme is continued until the condition on h is satisfied, or the Ap of the iteration 
becomes less than a given value. If the latter occurs, the iteration is terminated and an 
interpolation performed to obtain the desired thermodynamic quantities. 
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Iterative Procedure for Standing Normal Shock 

The conservation relations for mass, momentum, and energy for a standing normal 
shock are 


P 1 U 1 - P 2 U 2 

( 4 ) 

P 1 + P 1 U 1 2 = p 2 + P 2 U 2 2 

( 5 ) 

h, +~\J. 2 = ho +^U 9 2 
1 2 1 4 2 2 

( 6 ) 


Both the direct solution and inverse solution to the conservation relations are used herein. 

In the direct solution, the free-stream conditions appearing on the left side of equa- 
tions (4) to (6) are considered known. The can be determined from the energy rela- 
tion (eq. (6)) 

h ‘ = h i + l u 1 2 < 7 > 

If measured p t is not available, the relation (see ref. 1) 

p t = O^Uj 2 (8) 

is used to obtain an accurate estimate of p^.. From a known p^. or estimated value 
from equation (8) and known value of h^, a value of p t is obtained from FIND (2). The 
initial estimate of is taken to be 0.955p t . This value of Pg is used in equations (4) 
to (6) to obtain corresponding values of pg, hg, and I^. The P 2 and h 2 are used as 
inputs to FIND (2) to obtain a value of Pg. This Pg from FIND (2) is compared to the 
initial estimate of p£. If these values are not within 0.1 percent, the p ^ obtained from 
FIND (2) is used in equations (4) to (6) to obtain upgraded values of P 2 , hg, and U 2 . 

This method of successive approximations is repeated until successive values of P 2 
from FIND (2) are within 0.1 percent. 

The postshock stagnation conditions are determined by assuming that the flow 
region immediately behind the normal shock to the stagnation point is isentropic (that is, 

s t w o s 2 w o\ 

). Now, the variation in conditions from immediately downstream of the 

R R } 
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shock to the stagnation point is relatively small ^that is, 2 ~ y E t anc * ^2 * ^tj* 

Thus, this region may be considered to behave as a calorically perfect gas. The p t 
can be accurately estimated from the ideal-gas isentropic relation 


Pt = P2 




\ y E,2/ y E,2 _1 

M 2 2 ) 


( 9 ) 


This Pj. is used as an input to FEND (2), in conjunction with h^, to obtain a value of 
s t W 0 /R. If this s t W 0 ^R is not within 0.05 percent of s 2 W 0 /r, the p t is varied 
(upward or downward) 4 percent in 1 -percent increments. The values of s^y/ 0 jji cor- 
responding to each p t and the are obtained from FIND (2). The value of p t 
s t W G ® 2 ^o 

resulting in — - — = — - — is obtained by interpolation. 

R R 

In the inverse solution, the stagnation-point conditions are obtained from FIND (2) 
where inferred lq from stagnation-point heat-transfer rate measurement and measured 
P t serve as inputs. An estimate of p 2 is made and used with s t W Q /R as inputs to 
FIND (3) to obtain p 2 and h 2 . The corresponding U 2 is obtained from the energy 
relation (eq. (6)) 


U 2 =^2(ht-h 2 ) (10) 

In the procedures using the inverse solution, either p 1 or Pj is known. Thus, the 
free-stream quantities on the left side of equations (4) to (6) can be determined. The pj 
and Pj are used in FIND (1) to obtain a value of hj. This hj from FIND (1) is com-: 
pared to the calculated hj from equation (6) and, if not within 0.1 percent, the p 2 is 
varied and the procedure repeated. This numerical iteration on p 2 is continued until 
the convergence criterion on h^ is satisfied. 

Prediction of Stagnation- Point Heat-Transfer Rates 

Because measured q t is a basic input datum, it is necessary that a relatively 
accurate means for predicting q^. in C0 2 -N2-0 2 -Ar mixtures be available. Consider- 
able effort has been directed toward obtainment of expressions for predicting q t at high 
velocities in air. In comparison, the number of similar studies in arbitrary gas mixtures 
is relatively small. 

Recently, however, an analytical study of stagnation-point convective heat transfer 
to an axisymmetric blunt body was performed for a number of gas mixtures in thermo- 
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chemical equilibrium (ref. 6). This study demonstrated that the simple, empirical rela- 
tion of reference 7 for predicting q t in Ar-H2~N2-C02 gas mixtures yields results 
which are in good agreement with analytical results, provided the Np r w of the individ^ 
ual gases are approximately equal. For the gases considered herein, Np r w varies 
from approximately 0.67 for Ar to approximately 0.7 for N2 and O2 and to approximately 
0.72 for CO2 at the T w of interest (T w from room temperature to 1500° k). Hence, 

according to reference 6, the empirical relation of reference 7 should provide relatively 
accurate predictions of q t for the present gas mixtures. This relation of reference 7, 
which was also employed in the air study of reference 1, is 



where 


V s W 0 



and the are given in the following table: 


( 11 ) 


( 12 ) 


Species i 

^i 

co 2 

4.3102 x 10- 4 

n 2 

3.6285 

02 

4,8945 

Ar 

5.4788 


(The value of 17^ for O2 is not presented in ref. 7. This value was deduced by compari- 
son of the of table I in ref. 7 for air to the Kj for an air mixture of CO2, N2, O2, 
and Ar. The air-mixture composition was taken as 78.084 percent N2, 20.946 percent 
O2, 0.937 percent Ar, and 0.033 percent CO2 by volume (ref. 1).) The h w is obtained 
directly from ROGO, where the required temperature input to ROGO is T w . For the 
range of T w considered herein (300° K < T w < 1500° K), h w may be considered to be 
a function of T w only and a precise knowledge of p w is not necessary for obtaining 
h w accurately. For procedures where (p^ is an input or p t can be obtained from 

equation (8), the p w is equal to ^p^ or (p t ) . (fot the single procedure where 

is not an input, as well as and (Ui) m (that is, p^. cannot be obtained 

from eq. (8)), p^ is simply set equal to 1.01325 x 10 5 N/m2.^ 
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Procedures for Determining Free-Stream and Postshock Conditions 

The procedures for determining free-stream and post-normal-shock flow conditions 
are identified in the computer program by ITEST. For convenience, this method of iden- 
tification is employed in the following discussion. The basic measured inputs and itera- 
tive procedure for crossing the normal shock are given in the following table: 


ITEST 

Measured inputs 

Shock- crossing procedure 

1 

Pp u l> Pt 

Direct 

2 

Pp Pp ^t 

Inverse 

3 

Pp Pp 

Direct 

4 

Pi’ ^P Pt 

Direct 

5 

p v P t , q t 

Inverse 

6 . 

Pp p p ^t 

Direct 

7 

Pp Up q t 

Direct 


Of the 10 possible combinations of the input parameters pp Pp Up p t , and q t , three 
have been omitted from the present study. The procedure involving inputs Pp Up and 
p t was omitted because of the existence of a relatively accurate expression (eq. (8)) 
relating all three quantities. Omission of the other two procedures (up p t , q t , and pp 
Ui, q t ) was based on the error analysis performed in reference 1, which showed these 
procedures to be less acceptable than others. 

The individual procedures are now discussed briefly. 

ITEST = 1 .- A value of is determined from the measured Uj and p^. by 
using equation (8). The measured pj^ and this Pj are used as inputs to FIND (1) to 
obtain the corresponding free-stream thermodynamic properties. The direct iterative 
shock-crossing procedure is performed. The calculated p^. is compared to the mea- 
sured p t and if not within a prescribed tolerance (TOLPT), p^ is adjusted by using 
the relation 


h = (?i) 



prev (p t ) 


c,prev 


(13) 


This new value of p^ is used, in conjunction with (p^j , as input into FIND (1). 

The corresponding value of h^ obtained from FIND (1) is used in the direct shock- 
crossing solution. The ^p^ is again compared to . This procedure of upgrading 
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Pi according to equation (13) is continued until (Pt) c is within TOLPT of (Pt) m * (Note 

that Pj from eq. (8) was used only to obtain a relatively accurate first estimate of hj 
and is not involved in the final phase of upgrading P^-j 

ITEST =2 .- Since p^. and q^ are measured quantities and rg and T w are 
considered known, the is determined from equation (11). This and the mea- 
sured p t are used as inputs to FIND (2) to obtain the postshock stagnation conditions. 

The inverse iterative shock-crossing procedure is performed, where (Pi) m is the known 

free-stream flow quantity, to obtain the corresponding static conditions immediately 
behind the normal shock and the free-stream conditions. 


ITEST = 3 .- The measured p^ and p^ are used as inputs to FIND (1) to obtain 
the corresponding free-stream conditions. After the free-stream thermodynamic prop- 
erties are determined and Uj known, the postshock conditions are obtained from the 
direct shock-crossing procedure. 

ITEST =4 .- The Uj is determined from equation (8) and the' free-stream thermo- 
dynamic properties are obtained from FIND (1) with the measured pj and p^ as inputs. 
The is obtained from the direct shock- crossing procedure and compared to (Pt) m - 

If not within the desired tolerance (TOLPT), the Uj is upgraded according to 



This new value of Uj is used in the direct shock- crossing solution. The is 
upgraded according to equation (14) until (Pt) c is within TOLPT of (Pt) m * 

ITEST = 5 .- The procedure of ITEST = 5 is similar to that of ITEST = 2, the only 
difference being that p^ (instead of p^ is used as the required known free-stream flow 
quantity in the inverse shock-crossing solution. 

ITEST = 6 .- The free-stream thermodynamic quantities are obtained from FIND (1) 
with inputs (Pi) m and ^ Uj is obtained by combining equations (7), (8), 

and (11) to give 


+ 2(h! 




(15) 


The Ui from equation (15) is relatively accurate, the only uncertainty being in the value 
of 0.97 appearing in the last term. With a value of Uj, the postshock conditions are 
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determined as discussed in ITEST = 3. The value of 0.97 appearing in equation (15) is 
compared with 



(16) 


If not within 0.1 percent, the value 0.97 in equation (15) is replaced by Cj of equa- 
tion (16) and the calculation procedure repeated. 

ITEST = 7 .- In order to estimate a value of hj, equation (15) is rearranged to give 


hi = h w + 




(17) 


Again, the value of 0.97 appearing under the square root is the only uncertainty in this 
expression. The free-stream thermodynamic quantities are obtained from FIND (4) with 
and hj as inputs. The postshock conditions are determined as discussed in 


ITEST = 3. The value of 0.97 appearing in equation (17) is compared with Cj of equa- 
tion (16). If not within 0.1 percent, the value 0.97 in equation (17) is replaced by Cj of 
equation (16) and the procedure repeated. 


In many flow studies, the parameters free-stream Mach number and free-stream 
unit Reynolds number are of interest. These two parameters are determined in each 
procedure. The is obtained simply by dividing Uj by a*. The N Re ± is 

defined as hence must be determined. The viscosity of a gas mixture of 

n- components may be approximated by the relation (see ref. 8) 


where 



i + t*„ a 

k 1J Xi 




*. . 
i] 




+ 


Wi\l/ 2 

WjJ 


(18) 


(19) 
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In the present study, the number of components is taken to be four. These components 
correspond to the initial composition of Ar, CO 2 , N 2 , and O 2 [that is, the Wj and Wj 
of eq. (19) correspond to the W Q of the four gases under consideration^. 

The viscosity for Ar was obtained by applying a curve fit to the results of refer- 
ence 9, which includes data from a number of sources. The resulting expression is 

H = 3.33 x io- 7 T 0 - 739 (20) 


This expression yields values of fj. within 10 percent of the results illustrated in ref- 
erence 9, for the temperature range of 150° K to 5000° K. 

Simple expressions for the viscosity, in the form 


M = (a 0 


+ a-,T + a 0 T z x 10‘ 


( 21 ) 


were obtained for CO 2 , N 2 , and C >2 for temperature ranges of 150° K to 1000° K and 
1000° K to 5000° K. The constants [ag, a^, and a. 2 ) of equation (21) are presented in 
the following table, along with the source used in obtaining these constants: 


Gas 

AT 

a 0 

a l 

“2 

Ref. 

C0 2 

150 to 1000 

-9.7303 x 10" 1 

5.5222 x 10-1 

-1.6139 X 10-4 

10 

co 2 

1000 to 5000 

1.5029 x 102 

2.4890 x 10- 1 

-6.1747 X 10-® 

10, 11 

n 2 

150 to 1000 

2.2992 x 10 1 

5.5586 x 10- 1 

-1.8436 X lO" 4 

10, 12 

N 2 

1000 to 5000 

1.6388 x 10 2 

2.4601 x lO’ 1 

-7.7388 X lO” 6 

10, 13 

02 

150 to 1000 

1.9939 x 10 1 

6.6236 x 10' 1 

-2.0410 X lO" 4 

10 

°2 

1000 to 5000 

2.3030 x 10 2 

2.5597 x 10' 1 

-2.2643 x lO’ 6 

10, 14 


Viscosities for temperatures from 150° K to 5000° K were calculated for air 
( x Ar = 9 - 37 x 10 “ 3 > ^02 = 3.3x10-4, X Ng = 7.8084 x 10" 1 , and Xq 2 = 2.0946 x 10" 
by using equation (18) and compared to values from reference 1. This comparison showed 
that the values of fj. from equation (18) were within 3 percent of those of reference 1 for 
this temperature range. For mixtures where Ar or N 2 predominate, equation (18) should 
provide reasonably accurate values of fi for temperatures up to 5000° K. For mixtures 
where CO 2 or O 2 predominate, the fi of equation (18) will contain somewhat larger 
uncertainties at T > 2000° K because of dissociation of the CO 2 or O 2 . For example', 
reference 11 shows that for T > 1700° K, the MCO 2 ^ as a P ressure dependence because 
of dissociation of CO 2 . However, variation in pressure from 10 2 to 10® N/m2, for a given 
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temperature, results in a maximum uncertainty in McC>2 * ess than percent for 
temperatures less than 5000° K (ref. 11). Thus, equation (18) should provide values of 
[X, for mixtures where CO 2 or O 2 predominate, of sufficient accuracy for many flow 
studies. ^Eq. (18) was employed for the Mars model of ref. 11 (X^r = 0.32, X^c^ = 0.43, 
and X N g = 0.25) and a temperature range of 1000° K to 5000° K. These \i from 
eq. (18) were within 6 percent of those presented in ref. 11 for p = 10^ N/m^ for this 
temperature range .j 

DISCUSSION 

Because of the greater simplicity associated with procedure ITEST = 3, this pro- 
cedure was the first to be examined. A check case was run with air (78.084 percent N 2 , 
20.946 percent O 2 , 0.937 percent Ar, and 0.033 percent CO 2 by volume) as the test gas. 

The inputs pj, p^, and were the same values as used in obtaining the sample print- 
out of reference 1. Comparison of these results for ITEST = 3 to the sample printout 
case of reference 1 showed excellent agreement. The 15.24-km/sec entry trajectory 
cases presented in reference 1 were run with the present ITEST = 3 procedure and the 
preceding air model. For this trajectory, which encompassed a range of T^ from 
1500° K to 12 000° K, excellent agreement was observed between the present results and 
those of reference 1. Following the successful check of ITEST = 3, the remaining six 
procedures were run with a common check case. These check cases used the same free- 
stream inputs as ITEST = 3, but used the calculated postshock outputs p|. and q^ of 
ITEST = 3 as inputs. For these check cases, TOLPT = 0.001 and TOLQT = 0.005. 

The computed flow quantities for all six procedures were observed to be in agreement 
with those of ITEST = 3. 

Additional cases were run with ITEST = 3 for mixtures of 95 percent CO 2 and 
5 percent N 2 (representative of Venus) and 85 percent CO 2 , 1 percent N 2 , and 14 percent 
Ar (representative of Mars). The free -stream inputs were the same as those employed in 
the check case for ITEST = 3 and are representative of an expansion tube test in air with 
heated helium driver. The results of these additional cases are presented in the sample 
printouts of appendix C. 

In order to examine the relative total computer time (computational and peripheral) 
required by the various procedures, each procedure was run on a Control Data 6600 com- 
puter system. The representative expansion- tube test case for air used in the computer 
time study of reference 1 was also used in the present time study. The limits on temper- 
ature for FIND (I) (see appendix B) were as follows: 

TMAX1 = 1500° K 

TMAX2 = 8000° K 
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TMIN1 = 200° K 
TMIN2 = 5000° K 

The total computer times for each procedure were obtained for the case where all 26 spe- 
cies (see appendix B) were considered and for the case where the air model was simpli- 
fied to 10 species (e~, At, N, N + , N2, O, 0 + , O2, NO, NO+). The times for these two 
cases are given in the following table: 


ITEST 

Total computer time, sec 

26-species model 

10-species model 

1 

125 

70 

2 

250 

120 

3 

80 

45 

4 

110 

65 

5 

300 

125 

6 

80 

45 

7 

125 

70 


The reduction in computer time for the simplified air model is illustrated by this table. 
For the test case used in this time study, elimination of the 16 species between the two 
air models did not introduce significant variations (greater than 0.5 percent) in calculated 
free-stream and postshock conditions. (See printout of appendix C for 26-species air 
model and 10-species air model.) This result is not surprising since the percent by vol- 
ume of CO2 in the initial air composition for the 26- species model was relatively small 
(0.033 percent); hence elimination of CO2 from the initial composition should not signifi- 
cantly influence the calculated flow quantities. (Elimination of CO2 from the initial air 
composition corresponds to a reduction of eight species, since the C-element is elimi- 
nated.) Examination of the mole fractions at the stagnation-point conditions for the 
26 -species air model shows that the dependence on doubly ionized species along with O", 
Ar + , N2 + , C>2 + , and O2" is much less than the remaining species; thus, elimination of 
these species should also not significantly influence the calculated flow quantities for the 
test case being considered. Therefore, judicious selection of species can result in a 
substantial reduction in computer time without appreciable loss in accuracy. (See appen- 
dix B for discussion of variation in the number of species as required in the present com- 
puter program.) 
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CONCLUDING REMARKS 


Data reduction procedures for determining free-stream and post-normal-shock 
flow conditions are presented. These procedures are applicable to flows of CO 2 , N 2 , O 2 , 
Ar, or mixtures of these gases and include the effects of dissociation and ionization. The 
assumption of thermochemical equilibrium flow is made. Although derived primarily to 
meet the immediate needs of the expansion tube of the Langley hot gas radiation research 
facility, these procedures are applicable to any supersonic or hypersonic test facility 
using these gases or mixtures thereof. 

These procedures are based on the following flow parameters measured in the 
immediate vicinity of the test section: 

(1) Stagnation pressure behind normal shock 

(2) Free-stream static pressure 

(3) Stagnation-point heat-transfer rate 

(4) Free-stream velocity 

(5) Free-stream density 

The various combinations of measured flow parameters are identified herein as ITEST 
and are: 


ITEST Measured flow parameter 

1 (1), (2), (4) 

2 (1), (2), (3) 

3 (2), (4), (5) 

4 (1), (2), (5) 

5 (1), (3), (5) 

6 (2), (3), (5) 

7 (3), (4), (5) 

These seven procedures are incorporated into a single computer program written in 
FORTRAN IV language. 

Cases run with procedure ITEST = 3 for an air mixture were observed to be in 
good agreement with the results of the real-air study of NASA TN D-6618 for stagnation 
temperatures from 1500° K to 12 000° K. 
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A computer time study for each procedure showed that a substantial reduction in 
computer time, without appreciable loss in accuracy, may be realized in many cases by 
judicious selection of the number of species. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 2, 1972. 
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APPENDIX A 


LISTING OF COMPUTER PROGRAM 

The data reduction procedures for determining free-stream and post-normal-shock 
flow conditions for CO 2 , N 2 , O 2 , Ar, or mixtures of these gases are incorporated into a 
single computer program. This program is written in FORTRAN IV language for the 
Control Data 6600 computer system. Machine requirements are 70 000-octal locations 
of core storage. A listing of the main program (MILLER) and subroutines FIND (I) and 
ROGO, with comments, is reproduced in the following pages. 
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APPENDIX B 


COMPUTER PROGRAM INPUTS 


The inputs necessary to utilize the included computer program are presented in the 
following table: 


Program (FORTRAN) symbol 

Description 

ITEST 

Data reduction procedure to be used 

RUN 

Facility test number 

T0LPT 

Convergence criteria for iteration involving pj. 

T0LQT 

Convergence criteria for iteration involving 

TWE 

T W 

RN 

r g 

PC 02 

Measured partial pressure of CO 2 

PN2 

Measured partial pressure of N 2 

P02 

Measured partial pressure of O 2 

PAR 

Measured partial pressure of Ar 

PMIX 

Measured pressure of mixture 

P1M 

( p l)m 

PTM 

(Pt)m 

QTM 

(^t)m 

RH01M 

(Pl)m 

U1M 

(^l)m 

TMAX1 

Maximum Tj expected 

TMAX2 

Maximum T^ expected 

TMIN1 

Minimum Tj expected 

TMIN2 

Minimum T 2 expected 

NUMSP 

Number of species considered 

JINDX 

Number of components considered 

IAR 

Position in species array of Ar 

IN2 

Position in species array of N 2 

102 

Position in species array of O 2 

IC02 

Position in species array of CO 2 
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APPENDIX B - Continued 


The FORTRAN NAMELIST capability is used for the input data with INP as the NAMELIST 
name. A sample listing of the input cards for the 26-species air test case (ITEST = 3) is 
as follows: 

* INP I T£ST = 3 * RUN = 1 , P 1 M = 1 0 40 «UTm = 6' 1 00 * RHoTM=Yb0485 ♦ RN = . 0 1 27 , PM I X _ =Tr 
PC02 = • 00033 ' PN2= . 78084 .PO? =720 946 ^^=Td0937T^UMSP=26V J I NDX=5 * 1 AR = 2, " 
IN2=8« IQg= 14 « 1C02=26 S ■ ~ _ ~ ” 

The units for the inputs which are physical quantities are given in the section 
entitled "SYMBOLS." In using a specific ITEST, only those three measured inputs (com- 
binations of P1M, U1M, RH01M, PTM, and QTM) associated with this ITEST need be 
included. (The required measured inputs corresponding to a specific ITEST are given 
in the section entitled "Procedures for Determining Free-Stream and Postshock Condi- 
tions.") It should be noted that the condition 

PMIX = PAR + PN2 + P02 + PC02 

must be satisfied. If Ar, N 2 , O 2 , or CO 2 is not included in the test gas, then the partial 
pressure (PAR, PN2, P02, and PC02, respectively) need not be included as input. The 
IAR, IN2, 102, and IC02 need be included as input only if Ar, N 2 , O 2 , or CO 2 , respectively, 
is included in the species deck. (The computer cards containing the thermodynamic data 
for the various species are denoted as the species deck and are discussed subsequently.) 
The following table presents assigned values given to input quantities when these quantities 
are omitted from INP: 


Input quantity omitted 
(program symbol) 

Assigned value 

T0LPT 

0.001 

T0LQT 

0.005 

TWE 

300 

RN 

0.01 

TMAX1 

8 000 

TMAX2 

16 000 

TMIN1 

0 

TMIN2 

0 
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APPENDIX B - Continued 

The 26 species and five components considered are as follows: 
Species: 


e“ 

0 

C 

Ar 

o+ 

C+ 

Ar+ 

0++ 

C++ 

Ar ++ 

0" 

C' 

N 

°2 

CO 

N + 

°2 + 

C0+ 

N++ 

°2 _ 

CN 

n 2 

NO 

co 2 

n 2 + 

+ 

O 
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Components: 

e~ 

Ar 

N 

O 

C 

The thermodynamic data for these species, which are tabulated in reference 2, are 
read into the computer program from cards. A listing of these cards for all 26 species 
is as follows: 


44 


APPENDIX B — Continued 


El- i 0 0 54847000-03 00000000 + 00 00000000 + 00 

20000000+01 OOOOOOOO+nO 

A 30 0 O 39944000+02 00000000+00 00000000+00 

10000000+01 00000000+00 50000000+01 93144000+05 30000000+01 

93751000+05 lOOOOOCO+Ol 94554000+05 30000000+01 95400000+05 

30000000+01 10410200+06 12000000+02 10550000+06 80000000+01 

10615000+06 10000000+01 10705400+06 80000000+01 10722000+06 

40000000+01 10800000+06 40000000+01 11128000+06 16000000+02 

11175000+06 80000000+01 11290000+06 20000000+02 11355000+06 

16000000+02 11475000+^6 80000000+01 11500000+06 30000000+01 

11666000+06 12000000+02 11696000+06 80000000+01 11717000+06 

10000000+01 11756300+06 16000000+02 11853000+06 36000000+02 

11930000+06 56000000+02 12025000+06 20000000+02 12070000+06 

32000000+0 2 12175000+06 36000000+02 12220000+06 56000000+02 

12270000+06 20000000+02 12350000+06 28000000+02 12413600+06 

A+ 30 0 0 39943000+02 15202350+14 00000000+00 

60000000+01 00000000+00 20080000+01 10872300+06 20000000+02 

13240000+06 12000000+02 13480000+06 60000000+01 13860000+06 

28000000+02 14270000+06 60000000+01 14520000+06 12000000+02 

14765000+06 10000000+02 14875000+06 10000000+02 15000000+06 

12000000+02 15516000+06 30000000+02 15830000+06 12000000+02 

16100000+06 20000000+01 16730900+06 14000000+02 17060000+06 

26000000+02 17300000+06 38000000+02 17480000+06 60000000+01 

17970000+06 40000000+02 18300000+06 60000000+02 18600000+06 

62000000+02 19020000+06 34000000+02 19220000+06 12000000+03 

19400000+06 84000000+02 19600000+06 48000000+02 20000000+06 

18000000+02 20500000+06 60000000+01 20859300+06 22000000+03 

21000000+06 50000000+02 21500000+06 70000000+02 22000000+06 

A++ 30 0 0 39943000+02 41851700+14 00000000+00 

50000000+01 00000000+00 30000000+01 11120000+04 10000000+01 

15700000+04 50000000+01 14010000+05 10000000+01 33267000+05 

90000000+01 11440000+06 10000000+02 12800000+06 28000000+02 

14465000+06 15000000+02 15695000+06 50000000+01 17437500+06 

15000000+02 18200000+06 75000000+02 18950000+06 15000000+02 

19660000+06 50000000+01 20000000+06 18000000+02 20470000+06 

21000000+02 20830000+06 27000000+02 21080000+06 33000000+02 

21450000+06 46000000+02 22450000+06 90000000+01 23150000+06 

50000000+01 23500000+06 26000000+02 24000000+06 50000000+02 

24603600+06 91000000+02 25100000+06 24000000+02 25800000+06 

12700000+03 27000000+06 12100000+03 27850000+06 11600000+03 

28200000+06 11100000+03 28600000+06 52400000+03 30200000+06 

N 30 0 0 14008000+02 47072900+13 00000000+00 

40000000+01 00000000+00 10000000+02 19225000+05 60000000+01 

28840000+05 12000000+02 83330000+05 60000000+01 861S0000+05 
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12000000+02 88140000+05 20000000+01 93582000+05 20000000+02 

94800000+05 12000000+02 95500000+05 40000000+01 96752000+05 

10000000+02 96810000+05 60000000+01 97800000+05 10000000+02 

99660000+05 18000000+02 10400000+06 34000000+02 10470000+06 

56000000+02 10500000+06 22000000+02 10660000+06 16000000+02 

10720000+06 1 6000000+02 10760000+06 54000000+02 10960000+06 

12000000+02 10986000+06 10600000+03 11035000+06 13200000+03 

11160000+06 13600000+03 11220000+06 60000000+01 11231000+06 

68000000+02 11282000+06 30000000+02 11291000+06 20000000+01 

11360000+06 10000000+02 11390000+06 14000000+02 11700000+06 

N+ 30 0 0 1 4n07000+02 18726070+14 00000000+00 

90000000+01 00000000+00 50000000+01 15316000+05 10000000+01 

32687000+05 50000000+01 47168000+05 15000000+02 92245000+05 

90000000+01 10922000+06 50000000+01 14418900+06 12000000+02 

14900000+06 30000000+01 15513000+06 30000000+01 16461200+06 

18000000+02 16665000+06 30000000+01 16889300+06 90000000+01 

17062000+06 50000000+01 17421200+06 10000000+01 17827400+06 

41000000+02 18709000+06 16000000+02 18910000+06 30000000+01 

19012100+06 12000000+02 19720000+06 30000000+02 20300000+06 

21000000+02 20570000+06 15300000+03 21050000+06 12000000+02 

21482800+06 12000000+02 21800000+06 54600000+03 22050000+06 

37000000+02 22300000+06 15000000+02 22600000+06 15000000+02 

22800000+06 14000000+n2 23030000+06 50000000+01 23425000+06 

N++ 30 0 0 14007000+02 47288300+14 OOOOOOOO+OO 

60000000+01 COOOOOOO+OO 12000000+02 57280000+05 10000000+02 

10102600+06 20000000+01 13104400+06 60000000+01 14592000+06 

40000000+01 18680200+06 10000000+02 20307800+06 20000000+01 

22130200+06 60000000+01 23040700+06 60000000+01 24569000+06 

10000000+02 26724200+06 12000000+02 28765000+06 60000000+01 

29721000+06 20000000+01 30108800+06 32000000+02 31050000+06 

40000000+01 31422400+06 22000000+02 31755000+06 14000000+02 

32028800+06 10000000+02 32104000+06 30000000+02 32900000+06 

22000000+02 33330000+06 22000000+02 33629000+06 20000000+02 

33955000+06 48000000+02 34300000+06 24000000+02 34800000+06 

32000000+02 36400000+06 12000000+02 36862000+06 22000000+02 

374 0 000 O +06 52000000+02 37900000+06 20000000+01 38070000+06 

N2 6 1 0 28016000+02 00000000+00 20000000+01 

10000000+01 00000000+00 300000CD+0 1 49757000+05 60000000+01 

59314000+05 20000000+O1 68953000+05 13000000+01 70700000+05 

60000000+0 1 87984000+^5 

19980000+01 17000000-01 23580700+04 14190000+02 33 

14400000+01 1 3000000-nl 14603700+04 13890000+02 26 

16380000+01 18400000-01 17341100+04 14470000+02 29 

16182000+01 18300000-01 16937000+04 13830000+02 30 
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14800000+01 16400000-ni 15300000+04 12000000+02 31 

18258000+01 19700000-01 20351000+04 17080000+02 29 

N2+ 4 1 0 2801 5000+02 15033360+14 20000000+01 

2000000O+01 00000000+00 40000000+01 90200000+04 20000000+01 

25570000 + 05 20000000 + '"' 1 64550000 + 05 

19322000 + 01 20200000 — 01 22072300 + 04 16220000 + 02 32 

17220000+01 1 8000000— 0 1 19028400+04 14910000+02 31 

20830000+01 1 9500.000 — '"' 1 24198400 + 04 23190000 + 02 25 

16500000+01 50000000-01 20500000+04 14920000+02 30 

O 28 0 0 160 00000+02 24674100+13 00000000+00 

50000000+01 00000000+00 30000000+01 15900000+03 10000000+01 

22700000+03 50000000+01 15868000+05 10000000+01 33792000+05 

50000000+01 73768000+05 30000000+01 76795000+05 15000000+02 

86630000+05 90000000+01 88630000+05 50000000+01 95476000+05 

30000000+01 96226000+05 36000000+02 97420000+05 15000000+02 

97488000+05 15000000+02 99094000+05 90000000+01 99680000+05 

90000000+01 10000000+06 15000000+02 10114000+06 50000000+01 

10211600+06 30000000+01 10241200+06 50000000+01 10266^00+06 

25000000+02 10286500+06 15000000+02 10290800+06 90000C00+01 

10386900+06 15000000+02 10400000+06 25000000+02 10538500+06 

15000000+02 10540800+06 56000000+02 10600000+06 56000000+02 

1 0700000+06 

0+30 0 0 15999000+02 15603890+14 00000000+00 

40000000+01 00000000+00 10000000+02 26820000+05 60000000+01 

40465000+05 12000000+O2 12000000+06 10000000+02 16599000+06 

12000000+02 18540000+06 600Q00CD+01 18900000+06 20000000+01 

19571000+06 20000000+01 20394200+06 42000000+02 20760000+06 

26000000+02 21280000+06 20000000+01 22685100+06 24000000+02 

23000000+06 80000000+H2 23270000+06 16000000+02 23390000+06 

18000000+02 23960000+06 40000000+02 24550000+06 20000000+02 

24800000+06 60000000+01 25030000+06 32000000+02 25190000+06 

44000000+02 25400000+06 84000000+02 25550000+06 11600000+03 

25600000+06 18000000+02 25810000+06 10000000+02 25930000+06 

72000000+02 26150000+06 16600000+03 26550000+06 13400000+03 

27000000+06 14200000+O3 27650000+06 24000000+02 283000C0+06 

0++ 30 0 0 15999000+02 49479960+14 00000000+00 

90000000+01 00000000+00 50000000+01 20271000+05 10000000+01 

43184000+05 50000000+01 60312000+05 15000000+02 12005000+06 

90000000+01 14238300+06 50000000+01 18704900+06 30000000+01 

19708700+06 30000000+01 21045900+06 12000000+02 27000000+06 

90000000+01 28390000+06 26000000+02 29400000+06 14000000+02 

30300000+06 10000000+01 31380100+06 50000000+02 32700000+06 

10000000+02 33200000+06 15000000+02 33870000+06 10000000+01 

34330300+06 90000000+01 35020000+06 12000000+02 35750000+06 
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46000000+02 


36500000+06 

48000000+02 

37050000+06 

1 5800000+03 


38000000+06 


62000000+02 

39400000+06 

90000000+02 

39800000+06 


55800000+03 


40340000+06 

72000000+02 

42500000+06 

31 100000+03 


43000000+06 


48000000+02 

43800000+06 

1 5000000+0 2 

44271000+06 


0- 1 

0 

0 16001000+02 10541000+13 00000000+00 


60000000+0 1 


OOOOOOOO+oO 





02 6 

1 

0 32000000+02 00000000+00 20000000+01 I 


3000000O+0 1 


00000000+00 

20000000+01 

78824000+04 

1 0000000+0 1 


1 31 20900+05 


30000000+01 

357 1 3000+05 

1 0000000+01 

3621 2700+05 


30000000+0 1 


493631 00+05 





1 4457000 + 0 1 


1 5800000-01 

1 5803600+04 

12070000+02 

26 


14264000+0 1 


l 71 OOOOO-o 1 

15093000+04 

1 2900000+02 

29 


14004000+0 1 


1 8200000-O 1 

14326900+04 

1 3950000+02 

25 


1 0500000 + 0 1 


35700000-0 1 

8 1 900000+03 

22500000+02 

8 


82600000+00 


20500000-01 

65040000+03 

1 7030000+02 

9 


81900000+00 


1 1 000000-01 

70036000+03 

80000000+0 1 

21 


02+ 4 

1 

0 31999000+02 11628080+14 20000000+01 


40000000+0 1 


00000000+00 

80000000+0 1 

31500000+05 

40000000+0 1 


38300000+05 


40000000+01 

481 00000+05 




1 6722000+0 1 


19800000-01 

1 8764000+04 

1 6530000+02 

28 


1 1 047000 + 0 1 


15800000-^1 

1 0356900+04 

1 0390000+02 

24 


1 0617000+0 1 


19100000-01 

90000000+03 

1 3400000+02 

16 


12873000+01 


22100000-01 

1 1 967700+04 

1 7090000+02 

17 


02- 3 

1 

0 32001000+02 -96232000+12 20000000+01 


40000000+0 1 
24200000+05 


OOOOOOOO+OO 

40000000+0 1 

1 3400000+05 

40000000+0 1 


1 2000000 + 0 1 


1 6000000-0 1 

1 3000000+04 

1 4000000 + 02 

29 


97000000+00 


1 7000000-0 1 

99000000+03 

1 5000000+02 

16 


92000000+00 


27000000—n 1 

56000000+03 

1 3000000+02 

1 0 


NO 8 

1 

0 30O08000+02 89860000+12 10000000+01 


20000000+0 1 


00000000+00 

20000000 + 0 1 

1 2090000+03 

20000000+0 1 


44200000+05 


40000000+0 1 

45440000+05 

20000000+0 1 

53290000+05 


4000000O+0 1 


52376000+05 

20000000+0 1 

60860000+05 

40000000+0 1 


60020000+05 
1 7046000+0 1 


1 7800000-O1 

1 9040300+04 

1 3970000+02 

27 


1 7046000 + 0 1 


1 7800000-01 

1 9036800+04 

1 3970000 + 02 

33 


1 9972000+01 


1 9300000-0 1 

23748000+04 

1 6460000+02 

35 


1 1260000+01 


1 5200000-01 

1 0376800+04 

76030000+01 

33 


20026000+0 1 


2 1 800000-01 

23239000+04 

22885000+02 

25 


20020000+0 1 


30GOCOOO— 0 1 

23950000+04 

1 5C00000+02 

39 


19863000+01 


1 8200000-nl 

23736000+04 

1 5850000+02 

37 


1 3300000+0 1 


1 9000000-01 

1 21 66000+04 

1 5880000+02 

18 


N0+ 5 

1 

0 30O07000+02 98240300+13 10000000+01 


1 0000000+0 1 


00000000+00 

60000000+01 

39982000+05 

30000000+0 1 
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5852300O+05 60000000+01 72384000+05 20000000+01 73083800+05 

20020000+01 20200000-O 1 23771000+04 16350000+02 36 

16800000+01 1 9000000— o 1 17400000+04 14500000+02 29 

13300000+01 16000000-01 12200000+04 95000000+01 31 

12600000+01 17000000-01 11400000+04 76000000+01 37 

15870000+01 24000000-01 16089000+04 23300000+02 17 

C 3o 0 0 12011000+02 71123800+13 00000000+00 

90000000+01 00000000+00 50000000+01 10194000+05 10000000+01 

21648000+05 50000000+01 33735000+05 90000000+01 60360000+05 

30000000+01 61982000+05 15000000+02 64090000+05 30000000+01 

68858000+05 15000000+02 69700000+05 30000000+01 70744000+05 

90000000+01 71365000+05 50000000+01 72611000+05 10000000+01 

73976000+05 90000000+01 75256000+05 50000000+01 77681000+05 

90000000+01 78 1 30000+O5 21000000+02 78230000+05 18000000+02 

78320000+05 10000000+02 78600000+05 90000000+01 79318000+05 

18000000+02 80400000+05 12000000+02 81200000+05 50000000+01 

81770000+05 1 0000000+01 82252000+05 53000000+02 83800000+05 

10300000+03 84000000+05 30000000+02 84940000+05 60000000+01 

85400000+05 50300000+03 86400000+05 19000000+02 86500000+05 

C+ 30 0 0 12010000+02 17971820+14 00000000+00 

60000000+01 00000000+00 12000000+02 43030000+05 10000000+02 

74931000+05 20000000+01 96494000+05 60000000+01 10180000+06 

60000000+01 11065000+06 20000000+01 11490000+06 20000000+01 

11653800+06 14000000+02 11940000+06 60000000+01 13173100+06 

12000000+02 13580000+06 20000000+02 13600000+06 40000000+01 

14202400+06 10000000+02 14555100+06 10000000+02 15046500+06 

20000000+01 15723400+06 60000000+01 16252200+06 12000000+02 

16700000+06 10000000+02 16812400+06 20000000+02 16890000+06 

20000000+01 17334800+06 60000000+01 17529300+06 16000000+02 

17835000+06 32000000+02 17900000+06 26000000+02 18200000+06 

40000000+01 18468900+06 12000000+02 18645000+06 10000000+02 

18860000+06 30000000+02 19550000+06 20000000+02 19657200+06 

C++ 30 0 0 1 2o 1 0000+C2 41492250+14 00000000+00 

10000000+Cl 00000000+00 90000000+01 52360000+05 30000000+01 

10235100+06 90000000+01 13742000+06 50000000+01 14587500+06 

10000000+01 18252000+06 30000000+01 23816100+06 10000000+01 

24717000+06 30000000+01 25893100+06 90000000+01 25966200+06 

15000000+02 26996000+06 50000000+01 27684300+06 15000000+02 

30910000+06 10000000+01 31172100+06 12000000+02 31870000+06 

61000000+02 32255000+06 50000000+01 32421200+06 12000000+02 

32800000+06 31000000+02 33300000+06 15000000+02 33762600+06 

13000000+02 34000000+06 70000000+01 34136800+06 30000000+01 

34325600+06 30000000+02 34500000+06 60000000+02 34660000+06 

70000000+01 34800000+06 12000000+02 37660000+06 18000000+02 
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38 1 50000+06 


1 3000000+02 

38434500+06 

34000000+02 

38600000+06 

C- 1 

0 

0 1 2n 1 2000+02 58994400+13 00000000+00 

60000000+0 1 


OOOOOOOO+nO 




CO 7 

1 

0 28ol 1000+02 -1138 

1300+13 10000000+01 

10000000+0 1 


OOOOOOOO+oO 

60000000+01 

48473900+05 

30000000+0 1 

55380000+05 


60000000+0 1 

6 1 784600+05 

20000000+0 1 

64746500+05 

3000000^+0 1 


83831 000+05 

1 0000000+01 

86917800+05 


19313000+0 1 


1 7480000-01 

21 7021 00+04 

1 3461 000+02 

4 1 

1 6810000+0 1 


1 9300000-01 

1 7392500+04 

1 4470000 + 02 

29 

1 3310000+0 1 


1 6000000-0 1 

1 2 1 80000+04 

95000000+01 

31 

12615000+0 1 


1 7000000-01 

1 1 377900+04 

76240000+01 

37 

161 16000+01 


22290000-0 1 

151 56100+04 

1 7250500+02 

21 

20750000+0 1 


33000000-01 

2 1980000+04 

1 3476000+02 

40 

1961 0000+0 1 


27000000-n 1 

20820700+04 

12092000+02 

42 

C0+ 3 

1 

0 280 1 0000+02 12383670+14 10000000+01 

20000000+01 


00000000+00 

40000000+01 

20407500+05 

20000000+0 1 

45633500+05 
1 9772000+0 1 


1 8960000-01 

221 42400+04 

15164 000+02 

30 

1 5894000+0 1 


1 9420000-0 1 

1 5620600+04 

13530000+02 

28 

17999000+0 1 


30250000-O 1 

1 7341800+04 

27927000+02 

15 

CN 3 

1 

0 26019000+02 45605600+13 10000000+01 

20000000+0 1 
25797850+05 


00000000+00 

40000000+01 

91 145900+04 

20000000+01 

1 8996000+0 1 


1 7350000-01 

20687050+04 

1 3144000+02 

29 

17165000+0 1 


1 7460000-0 1 

1 8144300+04 

1 2883000+02 

34 

19701000+01 


22150000-01 

21641300+04 

20250000+02 

26 

C02 1 

2 

0 44nll000+02 -39314600+13 20000000+01 

1 0000000+0 1 


OOOOOOOO+oO 




39060000+00 


22000000-02 




66733000+03 


00000000+00 

66733000+03 

00000000+00 

1 3428600+04 

00000000+00 


23493000+04 

00000000+00 




Cases may arise in which the user of this program may be interested in a gas or 
gas mixture involving fewer than 26 species. If so, only the species cards for the species 
being considered need be included. The number of species (NUMSP) and number of com- 
ponents (JINDX) are basic inputs and must be adjusted accordingly. The inputs IAR, IN2, 
102, and IC02, which designate the position of Ar, N 2 , O 2 , and CO 2 , respectively, in the 
species array, must also be adjusted to provide meaningful initial estimates for the spe- 
cies concentrations as required by ROGO. Finally, the number of atoms of a particular 
component per particle of a species A(I,J) must be designated. For electrons, a charge 
constraint is used instead of a mass constraint. 

The modifications required when the number of species is varied is best illustrated 
by example. First, consider a C02-N2-C>2-Ar mixture in which double ionization is 
included (that is, all 26 species are included). A suggested aid to this case is first to 
construct a table of A(I,J) as follows: 
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The sequence 0, 1, 0, 0, 0 for species Ar denotes that the only component contained in 
this species is Ar. (A zero denotes that a particular component does not appear in the 
species being considered.) In the sequence -1, 1, 0, 0, 0 for Ar+, the -1 denotes that Ar+ 
is a positive ion, and so forth. 

Now, for the mixture being considered, NUMSP = 26, JINDX = 5, IAR = 2, 

IN2 = 8, 102 = 14, and IC02 = 26. Naturally, all 26 species decks are included and are 
placed in the order that they appear in the preceding table of A(I,J). It now remains to 
read in the A(I,J) as illustrated by the following listing: 

1 ooooo 1 000 -1 1 000 

-2 1 0 0 0 0 0 1 00-1 0 1 00 
-2 0 1 0000200 -1 0200 
000 10 -10010 -2 0010 

100 1 000020 -1 0020 

1002000110 -10110 
00001 -1 0001 -2 0001 

10001 00011-1 001 1 

0 0 1 0 1 0 0 0 2 1 

Thus, the first card contains the A(I,J) for e", Ar, and Ar + ; the second, for Ar ++ , N, and 
N + ; and so forth. 

The sequence of cards following the main program and subroutines is as follows: 
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If the simplified air model used in the previously discussed computer time study is 
now considered, the table of A(I,J) takes the following form: 


12 3 4 



^\Component (J) 
Species 

e" 

Ar 

N 

0 

1 

e" 

1 

0 

0 

0 

2 

Ar 

0 

1 

0 

0 

3 

N 

0 

0 

1 

0 

4 

N + 

-1 

0 

1 

0 

5 

n 2 

0 

0 

2 

0 

6 

O 

0 

0 

0 

1 

7 

o+ 

-1 

0 

0 

1 

8 

02 

0 

0 

0 

2 

9 

NO 

0 

0 

1 

1 

10 

NO+ 

-1 

0 

1 

1 


Then, NUMSP = 10, JINDX = 4, IAR = 2, IN2 = 5, and 102 = 8. Only the decks for the 
10 species considered are included, and these are placed in the order that they appear in 
the preceding table of A(I,J). The A(I,J) are read in as illustrated by the following listing: 


lOOOO 
0 0 0 2 0 
0 2 0 0 1 


1 0 0 0 0 1 0 -1 0 1 
000 1-1 00 1 00 
1-1011 


Since the format for reading in the A(I,J) is the same as for the 26-species case, the first 
card contains all the A(I,J) for e", Ar, N, and the first three A(I,J) for N + ; the second card 
contains the fourth A(I,J) for N+, all the A(I,J) for N 2 , O, and 0 + , and the first two A(I,J) 
for O 2 ; and so forth. 
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APPENDIX C 


SAMPLE DATA PRINTOUTS OF COMPUTER PROGRAM 

Sample data printouts are given herein. These printouts are for a 26 -species air 
model, a 10-species air model, a 95-percent CO 2 , 5-percent N 2 Venus model, and an 
85-percent CO 2 , 1-percent N 2 , 14-percent Ar Mars model, respectively. Procedure 
ITEST = 3 was used to obtain these printouts. 

The headings for the various free-stream and postshock flow conditions of interest 
correspond to those in the section entitled "SYMBOLS." The units of these flow condi- 
tions are also given in this section. The sample computer printouts are as follows: 
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